Epigenetic profiling reveals key genes and cis-regulatory networks specific to human parathyroids.

In all terrestrial vertebrates, the parathyroid glands are critical regulators of calcium homeostasis and the sole source of parathyroid hormone (PTH). Hyperparathyroidism and hypoparathyroidism are clinically important disorders affecting multiple organs. However, our knowledge regarding regulatory mechanisms governing the parathyroids has remained limited. Here, we present the comprehensive maps of the chromatin landscape of the human parathyroid glands, identifying active regulatory elements and chromatin interactions. These data allow us to define regulatory circuits and previously unidentified genes that play crucial roles in parathyroid biology. We experimentally validate candidate parathyroid-specific enhancers and demonstrate their integration with GWAS SNPs for parathyroid-related diseases and traits. For instance, we observe reduced activity of a parathyroid-specific enhancer of the Calcium Sensing Receptor gene, which contains a risk allele associated with higher PTH levels compared to the wildtype allele. Our datasets provide a valuable resource for unraveling the mechanisms governing parathyroid gland regulation in health and disease.

Activation of CD8+ T Cells in Chronic Obstructive Pulmonary Disease Lung.

Rationale: Despite the importance of inflammation in chronic obstructive pulmonary disease (COPD), the immune cell landscape in the lung tissue of patients with mild-moderate disease has not been well characterized at the single-cell and molecular level. Objectives: To define the immune cell landscape in lung tissue from patients with mild-moderate COPD at single-cell resolution. Methods: We performed single-cell transcriptomic, proteomic, and T-cell receptor repertoire analyses on lung tissue from patients with mild-moderate COPD (n = 5, Global Initiative for Chronic Obstructive Lung Disease I or II), emphysema without airflow obstruction (n = 5), end-stage COPD (n = 2), control (n = 6), or donors (n = 4). We validated in an independent patient cohort (N = 929) and integrated with the Hhip+/- murine model of COPD. Measurements and Main Results: Mild-moderate COPD lungs have increased abundance of two CD8+ T cell subpopulations: cytotoxic KLRG1+TIGIT+CX3CR1+ TEMRA (T effector memory CD45RA+) cells, and DNAM-1+CCR5+ T resident memory (TRM) cells. These CD8+ T cells interact with myeloid and alveolar type II cells via IFNG and have hyperexpanded T-cell receptor clonotypes. In an independent cohort, the CD8+KLRG1+ TEMRA cells are increased in mild-moderate COPD lung compared with control or end-stage COPD lung. Human CD8+KLRG1+ TEMRA cells are similar to CD8+ T cells driving inflammation in an aging-related murine model of COPD. Conclusions: CD8+ TEMRA cells are increased in mild-moderate COPD lung and may contribute to inflammation that precedes severe disease. Further study of these CD8+ T cells may have therapeutic implications for preventing severe COPD.

Cigarette smoking-associated isoform switching and 3' UTR lengthening via alternative polyadenylation.

Cigarette smoking induces a profound transcriptomic and systemic inflammatory response. Previous studies have focused on gene level differential expression of smoking, but the genome-wide effects of smoking on alternative isoform regulation have not yet been described. We conducted RNA sequencing in whole-blood samples of 454 current and 767 former smokers in the COPDGene Study, and we analyzed the effects of smoking on differential usage of isoforms and exons. At 10% FDR, we detected 3167 differentially expressed genes, 945 differentially used isoforms and 160 differentially used exons. Isoform switch analysis revealed widespread 3' UTR lengthening associated with cigarette smoking. The lengthening of these 3' UTRs was consistent with alternative usage of distal polyadenylation sites, and these extended 3' UTR regions were significantly enriched with functional sequence elements including microRNA and RNA-protein binding sites. These findings warrant further studies on alternative polyadenylation events as potential biomarkers and novel therapeutic targets for smoking-related diseases.

Essential histone chaperones collaborate to regulate transcription and chromatin integrity.

Histone chaperones are critical for controlling chromatin integrity during transcription, DNA replication, and DNA repair. Three conserved and essential chaperones, Spt6, Spn1/Iws1, and FACT, associate with elongating RNA polymerase II and interact with each other physically and/or functionally; however, there is little understanding of their individual functions or their relationships with each other. In this study, we selected for suppressors of a temperature-sensitive spt6 mutation that disrupts the Spt6-Spn1 physical interaction and that also causes both transcription and chromatin defects. This selection identified novel mutations in FACT. Surprisingly, suppression by FACT did not restore the Spt6-Spn1 interaction, based on coimmunoprecipitation, ChIP, and mass spectrometry experiments. Furthermore, suppression by FACT bypassed the complete loss of Spn1. Interestingly, the FACT suppressor mutations cluster along the FACT-nucleosome interface, suggesting that they alter FACT-nucleosome interactions. In agreement with this observation, we showed that the spt6 mutation that disrupts the Spt6-Spn1 interaction caused an elevated level of FACT association with chromatin, while the FACT suppressors reduced the level of FACT-chromatin association, thereby restoring a normal Spt6-FACT balance on chromatin. Taken together, these studies reveal previously unknown regulation between histone chaperones that is critical for their essential in vivo functions.

HiTea: a computational pipeline to identify non-reference transposable element insertions in Hi-C data.

Hi-C is a common technique for assessing 3D chromatin conformation. Recent studies have shown that long-range interaction information in Hi-C data can be used to generate chromosome-length genome assemblies and identify large-scale structural variations. Here, we demonstrate the use of Hi-C data in detecting mobile transposable element (TE) insertions genome-wide. Our pipeline Hi-C-based TE analyzer (HiTea) capitalizes on clipped Hi-C reads and is aided by a high proportion of discordant read pairs in Hi-C data to detect insertions of three major families of active human TEs. Despite the uneven genome coverage in Hi-C data, HiTea is competitive with the existing callers based on whole-genome sequencing (WGS) data and can supplement the WGS-based characterization of the TE-insertion landscape. We employ the pipeline to identify TE-insertions from human cell-line Hi-C samples. AVAILABILITY AND IMPLEMENTATION: HiTea is available at https://github.com/parklab/HiTea and as a Docker image. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The conserved elongation factor Spn1 is required for normal transcription, histone modifications, and splicing in Saccharomyces cerevisiae.

Spn1/Iws1 is a conserved protein involved in transcription and chromatin dynamics, yet its general in vivo requirement for these functions is unknown. Using a Spn1 depletion system in Saccharomyces cerevisiae, we demonstrate that Spn1 broadly influences several aspects of gene expression on a genome-wide scale. We show that Spn1 is globally required for normal mRNA levels and for normal splicing of ribosomal protein transcripts. Furthermore, Spn1 maintains the localization of H3K36 and H3K4 methylation across the genome and is required for normal histone levels at highly expressed genes. Finally, we show that the association of Spn1 with the transcription machinery is strongly dependent on its binding partner, Spt6, while the association of Spt6 and Set2 with transcribed regions is partially dependent on Spn1. Taken together, our results show that Spn1 affects multiple aspects of gene expression and provide additional evidence that it functions as a histone chaperone in vivo.

Epigenetic transcriptional reprogramming by WT1 mediates a repair response during podocyte injury.

In the context of human disease, the mechanisms whereby transcription factors reprogram gene expression in reparative responses to injury are not well understood. We have studied the mechanisms of transcriptional reprogramming in disease using murine kidney podocytes as a model for tissue injury. Podocytes are a crucial component of glomeruli, the filtration units of each nephron. Podocyte injury is the initial event in many processes that lead to end-stage kidney disease. Wilms tumor-1 (WT1) is a master regulator of gene expression in podocytes, binding nearly all genes known to be crucial for maintenance of the glomerular filtration barrier. Using murine models and human kidney organoids, we investigated WT1-mediated transcriptional reprogramming during the course of podocyte injury. Reprogramming the transcriptome involved highly dynamic changes in the binding of WT1 to target genes during a reparative injury response, affecting chromatin state and expression levels of target genes.

HiNT: a computational method for detecting copy number variations and translocations from Hi-C data.

The three-dimensional conformation of a genome can be profiled using Hi-C, a technique that combines chromatin conformation capture with high-throughput sequencing. However, structural variations often yield features that can be mistaken for chromosomal interactions. Here, we describe a computational method HiNT (Hi-C for copy Number variation and Translocation detection), which detects copy number variations and interchromosomal translocations within Hi-C data with breakpoints at single base-pair resolution. We demonstrate that HiNT outperforms existing methods on both simulated and real data. We also show that Hi-C can supplement whole-genome sequencing in structure variant detection by locating breakpoints in repetitive regions.

Comprehensive analysis of chromothripsis in 2,658 human cancers using whole-genome sequencing.

Chromothripsis is a mutational phenomenon characterized by massive, clustered genomic rearrangements that occurs in cancer and other diseases. Recent studies in selected cancer types have suggested that chromothripsis may be more common than initially inferred from low-resolution copy-number data. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we analyze patterns of chromothripsis across 2,658 tumors from 38 cancer types using whole-genome sequencing data. We find that chromothripsis events are pervasive across cancers, with a frequency of more than 50% in several cancer types. Whereas canonical chromothripsis profiles display oscillations between two copy-number states, a considerable fraction of events involve multiple chromosomes and additional structural alterations. In addition to non-homologous end joining, we detect signatures of replication-associated processes and templated insertions. Chromothripsis contributes to oncogene amplification and to inactivation of genes such as mismatch-repair-related genes. These findings show that chromothripsis is a major process that drives genome evolution in human cancer.

CHRAC/ACF contribute to the repressive ground state of chromatin.

The chromatin remodeling complexes chromatin accessibility complex and ATP-utilizing chromatin assembly and remodeling factor (ACF) combine the ATPase ISWI with the signature subunit ACF1. These enzymes catalyze well-studied nucleosome sliding reactions in vitro, but how their actions affect physiological gene expression remains unclear. Here, we explored the influence of Drosophila melanogaster chromatin accessibility complex/ACF on transcription by using complementary gain- and loss-of-function approaches. Targeting ACF1 to multiple reporter genes inserted at many different genomic locations revealed a context-dependent inactivation of poorly transcribed reporters in repressive chromatin. Accordingly, single-embryo transcriptome analysis of an Acf knock-out allele showed that only lowly expressed genes are derepressed in the absence of ACF1. Finally, the nucleosome arrays in Acf-deficient chromatin show loss of physiological regularity, particularly in transcriptionally inactive domains. Taken together, our results highlight that ACF1-containing remodeling factors contribute to the establishment of an inactive ground state of the genome through chromatin organization.

Genome-wide Rules of Nucleosome Phasing in Drosophila.

Regular successions of positioned nucleosomes, or phased nucleosome arrays (PNAs), are predominantly known from transcriptional start sites (TSSs). It is unclear whether PNAs occur elsewhere in the genome. To generate a comprehensive inventory of PNAs for Drosophila, we applied spectral analysis to nucleosome maps and identified thousands of PNAs throughout the genome. About half of them are not near TSSs and are strongly enriched for an uncharacterized sequence motif. Through genome-wide reconstitution of physiological chromatin in Drosophila embryo extracts, we uncovered the molecular basis of PNA formation. We identified Phaser, an unstudied zinc finger protein that positions nucleosomes flanking the motif. It also revealed how the global activity of the chromatin remodelers CHRAC/ACF, together with local barrier elements, generates islands of regular phasing throughout the genome. Our work demonstrates the potential of chromatin assembly by embryo extracts as a powerful tool to reconstitute chromatin features on a global scale in vitro.

EED, a member of the polycomb group, is required for nephron differentiation and the maintenance of nephron progenitor cells.

Epigenetic regulation of gene expression has a crucial role allowing for the self-renewal and differentiation of stem and progenitor populations during organogenesis. The mammalian kidney maintains a population of self-renewing stem cells that differentiate to give rise to thousands of nephrons, which are the functional units that carry out filtration to maintain physiological homeostasis. The polycomb repressive complex 2 (PRC2) epigenetically represses gene expression during development by placing the H3K27me3 mark on histone H3 at promoter and enhancer sites, resulting in gene silencing. To understand the role of PRC2 in nephron differentiation, we conditionally inactivated the Eed gene, which encodes a nonredundant component of the PRC2 complex, in nephron progenitor cells. Resultant kidneys were smaller and showed premature loss of progenitor cells. The progenitors in Eed mutant mice that were induced to differentiate did not develop into properly formed nephrons. Lhx1, normally expressed in the renal vesicle, was overexpressed in kidneys of Eed mutant mice. Thus, PRC2 has a crucial role in suppressing the expression of genes that maintain the progenitor state, allowing nephron differentiation to proceed.

A role for tuned levels of nucleosome remodeler subunit ACF1 during Drosophila oogenesis.

The Chromatin Accessibility Complex (CHRAC) consists of the ATPase ISWI, the large ACF1 subunit and a pair of small histone-like proteins, CHRAC-14/16. CHRAC is a prototypical nucleosome sliding factor that mobilizes nucleosomes to improve the regularity and integrity of the chromatin fiber. This may facilitate the formation of repressive chromatin. Expression of the signature subunit ACF1 is restricted during embryonic development, but remains high in primordial germ cells. Therefore, we explored roles for ACF1 during Drosophila oogenesis. ACF1 is expressed in somatic and germline cells, with notable enrichment in germline stem cells and oocytes. The asymmetrical localization of ACF1 to these cells depends on the transport of the Acf1 mRNA by the Bicaudal-D/Egalitarian complex. Loss of ACF1 function in the novel Acf1(7) allele leads to defective egg chambers and their elimination through apoptosis. In addition, we find a variety of unusual 16-cell cyst packaging phenotypes in the previously known Acf1(1) allele, with a striking prevalence of egg chambers with two functional oocytes at opposite poles. Surprisingly, we found that the Acf1(1) deletion--despite disruption of the Acf1 reading frame--expresses low levels of a PHD-bromodomain module from the C-terminus of ACF1 that becomes enriched in oocytes. Expression of this module from the Acf1 genomic locus leads to packaging defects in the absence of functional ACF1, suggesting competitive interactions with unknown target molecules. Remarkably, a two-fold overexpression of CHRAC (ACF1 and CHRAC-16) leads to increased apoptosis and packaging defects. Evidently, finely tuned CHRAC levels are required for proper oogenesis.

Active promoters give rise to false positive 'Phantom Peaks' in ChIP-seq experiments.

Chromatin immunoprecipitation (ChIP) is widely used to identify chromosomal binding sites. Chromatin proteins are cross-linked to their target sequences in living cells. The purified chromatin is sheared and the relevant protein is enriched by immunoprecipitation with specific antibodies. The co-purifying genomic DNA is then determined by massive parallel sequencing (ChIP-seq).We applied ChIP-seq to map the chromosomal binding sites for two ISWI-containing nucleosome remodeling factors, ACF and RSF, in Drosophila embryos. Employing several polyclonal and monoclonal antibodies directed against their signature subunits, ACF1 and RSF-1, robust profiles were obtained indicating that both remodelers co-occupied a large set of active promoters.Further validation included controls using chromatin of mutant embryos that do not express ACF1 or RSF-1. Surprisingly, the ChIP-seq profiles were unchanged, suggesting that they were not due to specific immunoprecipitation. Conservative analysis lists about 3000 chromosomal loci, mostly active promoters that are prone to non-specific enrichment in ChIP and appear as 'Phantom Peaks'. These peaks are not obtained with pre-immune serum and are not prominent in input chromatin.Mining the modENCODE ChIP-seq profiles identifies potential Phantom Peaks in many profiles of epigenetic regulators. These profiles and other ChIP-seq data featuring prominent Phantom Peaks must be validated with chromatin from cells in which the protein of interest has been depleted.

Molecular analysis of gut microbiota in obesity among Indian individuals.

Obesity is a consequence of a complex interplay between the host genome and the prevalent obesogenic factors among the modern communities. The role of gut microbiota in the pathogenesis of the disorder was recently discovered; however, 16S-rRNA-based surveys revealed compelling but community-specific data. Considering this, despite unique diets, dietary habits and an uprising trend in obesity, the Indian counterparts are poorly studied. Here, we report a comparative analysis and quantification of dominant gut microbiota of lean, normal, obese and surgically treated obese individuals of Indian origin. Representative gut microbial diversity was assessed by sequencing fecal 16S rRNA libraries for each group (n=5) with a total of over 3000 sequences. We detected no evident trend in the distribution of the predominant bacterial phyla, Bacteroidetes and Firmicutes. At the genus level, the bacteria of genus Bacteroides were prominent among the obese individuals, which was further confirmed by qPCR (P less than 0.05). In addition, a remarkably high archaeal density with elevated fecal SCFA levels was also noted in the obese group. On the contrary, the treated-obese individuals exhibited comparatively reduced Bacteroides and archaeal counts along with reduced fecal SCFAs. In conclusion, the study successfully identified a representative microbial diversity in the Indian subjects and demonstrated the prominence of certain bacterial groups in obese individuals; nevertheless, further studies are essential to understand their role in obesity.